SIRAH Late Harvest: Coarse-Grained Models for Protein Glycosylation.

Glycans constitute one of the most complex families of biological molecules. Despite their crucial role in a plethora of biological processes, they remain largely uncharacterized because of their high complexity. Their intrinsic flexibility and the vast variability associated with the many combination possibilities have hampered their experimental determination. Although theoretical methods have proven to be a valid alternative to the study of glycans, the large size associated with polysaccharides, proteoglycans, and glycolipids poses significant challenges to a fully atomistic description of biologically relevant glycoconjugates. On the other hand, the exquisite dependence on hydrogen bonds to determine glycans' structure makes the development of simplified or coarse-grained (CG) representations extremely challenging. This is particularly the case when glycan representations are expected to be compatible with CG force fields that include several molecular types. We introduce a CG representation able to simulate a wide variety of polysaccharides and common glycosylation motifs in proteins, which is fully compatible with the CG SIRAH force field. Examples of application to N-glycosylated proteins, including antibody recognition and calcium-mediated glycan-protein interactions, highlight the versatility of the enlarged set of CG molecules provided by SIRAH.

AmberTools.

AmberTools is a free and open-source collection of programs used to set up, run, and analyze molecular simulations. The newer features contained within AmberTools23 are briefly described in this Application note.

New Insights into the Role of the Trypanosoma cruzi Aldo-Keto Reductase TcAKR.

Chagas disease is a zoonotic infectious disease caused by the protozoan parasite Trypanosoma cruzi. It is distributed worldwide, affecting around 7 million people; there is no effective treatment, and it constitutes a leading cause of disability and premature death in the Americas. Only two drugs are currently approved for the treatment, Benznidazole and Nifurtimox, and both have to be activated by reducing the nitro-group. The T. cruzi aldo-keto reductase (TcAKR) has been related to the metabolism of benznidazole. TcAKR has been extensively studied, being most efforts focused on characterizing its implication in trypanocidal drug metabolism; however, little is known regarding its biological role. Here, we found that TcAKR is confined, throughout the entire life cycle, into the parasite mitochondria providing new insights into its biological function. In particular, in epimastigotes, TcAKR is associated with the kinetoplast, which suggests additional roles of the protein. The upregulation of TcAKR, which does not affect TcOYE expression, was correlated with an increase in PGF2α, suggesting that this enzyme is related to PGF2α synthesis in T. cruzi. Structural analysis showed that TcAKR contains a catalytic tetrad conserved in the AKR superfamily. Finally, we found that TcAKR is also involved in Nfx metabolization.

Thiol-based chemical probes exhibit antiviral activity against SARS-CoV-2 via allosteric disulfide disruption in the spike glycoprotein.

The development of small-molecules targeting different components of SARS-CoV-2 is a key strategy to complement antibody-based treatments and vaccination campaigns in managing the COVID-19 pandemic. Here, we show that two thiol-based chemical probes that act as reducing agents, P2119 and P2165, inhibit infection by human coronaviruses, including SARS-CoV-2, and decrease the binding of spike glycoprotein to its receptor, the angiotensin-converting enzyme 2 (ACE2). Proteomics and reactive cysteine profiling link the antiviral activity to the reduction of key disulfides, specifically by disruption of the Cys379-Cys432 and Cys391-Cys525 pairs distal to the receptor binding motif in the receptor binding domain (RBD) of the spike glycoprotein. Computational analyses provide insight into conformation changes that occur when these disulfides break or form, consistent with an allosteric role, and indicate that P2119/P2165 target a conserved hydrophobic binding pocket in the RBD with the benzyl thiol-reducing moiety pointed directly toward Cys432. These collective findings establish the vulnerability of human coronaviruses to thiol-based chemical probes and lay the groundwork for developing compounds of this class, as a strategy to inhibit the SARS-CoV-2 infection by shifting the spike glycoprotein redox scaffold.

Emergence and Spread of a B.1.1.28-Derived P.6 Lineage with Q675H and Q677H Spike Mutations in Uruguay.

Uruguay controlled the viral dissemination during the first nine months of the SARS-CoV-2 pandemic. Unfortunately, towards the end of 2020, the number of daily new cases exponentially increased. Herein, we analyzed the country-wide genetic diversity of SARS-CoV-2 between November 2020 and April 2021. We identified that the most prevalent viral variant during the first epidemic wave in Uruguay (December 2020-February 2021) was a B.1.1.28 sublineage carrying Spike mutations Q675H + Q677H, now designated as P.6, followed by lineages P.2 and P.7. P.6 probably arose around November 2020, in Montevideo, Uruguay's capital department, and rapidly spread to other departments, with evidence of further local transmission clusters; it also spread sporadically to the USA and Spain. The more efficient dissemination of lineage P.6 with respect to P.2 and P.7 and the presence of mutations (Q675H and Q677H) in the proximity of the key cleavage site at the S1/S2 boundary suggest that P.6 may be more transmissible than other lineages co-circulating in Uruguay. Although P.6 was replaced by the variant of concern (VOC) P.1 as the predominant lineage in Uruguay since April 2021, the monitoring of the concurrent emergence of Q675H + Q677H in VOCs should be of worldwide interest.

Fighting viruses with computers, right now.

The synergistic conjunction of various technological revolutions with the accumulated knowledge and workflows is rapidly transforming several scientific fields. Particularly, Virology can now feed from accurate physical models, polished computational tools, and massive computational power to readily integrate high-resolution structures into biological representations of unprecedented detail. That preparedness allows for the first time to get crucial information for vaccine and drug design from in-silico experiments against emerging pathogens of worldwide concern at relevant action windows. The present work reviews some of the main milestones leading to these breakthroughs in Computational Virology, providing an outlook for future developments in capacity building and accessibility to computational resources.

CUTie2: The Attack of the Cyclic Nucleotide Sensor Clones.

The detection of small molecules in living cells using genetically encoded FRET sensors has revolutionized our understanding of signaling pathways at the sub-cellular level. However, engineering fluorescent proteins and specific binding domains to create new sensors remains challenging because of the difficulties associated with the large size of the polypeptides involved, and their intrinsically huge conformational variability. Indeed, FRET sensors' design still relies on vague structural notions, and trial and error combinations of linkers and protein modules. We recently designed a FRET sensor for the second messenger cAMP named CUTie (Cyclic nucleotide Universal Tag for imaging experiments), which granted sub-micrometer resolution in living cells. Here we apply a combination of sequence/structure analysis to produce a new-generation FRET sensor for the second messenger cGMP based on Protein kinase G I (PKGI), which we named CUTie2. Coarse-grained molecular dynamics simulations achieved an exhaustive sampling of the relevant spatio-temporal coordinates providing a quasi-quantitative prediction of the FRET efficiency, as confirmed by in vitro experiments. Moreover, biochemical characterization showed that the cGMP binding module maintains virtually the same affinity and selectivity for its ligand thant the full-length protein. The computational approach proposed here is easily generalizable to other allosteric protein modules, providing a cost effective-strategy for the custom design of FRET sensors.

Wrapping Up Viruses at Multiscale Resolution: Optimizing PACKMOL and SIRAH Execution for Simulating the Zika Virus.

Simulating huge biomolecular complexes of million atoms at relevant biological time scales is becoming accessible to the broad scientific community. That proves to be crucial for urgent responses against emergent diseases in real time. Yet, there are still issues to sort regarding the system setup so that molecular dynamics (MD) simulations can be run in a simple and standard way. Here, we introduce an optimized pipeline for building and simulating enveloped virus-like particles (VLP). First, the membrane packing problem is tackled with new features and optimized options in PACKMOL. This allows preparing accurate membrane models of thousands of lipids in the context of a VLP within a few hours using a single CPU. Then, the assembly of the VLP system is done within the multiscale framework of the coarse-grained SIRAH force field. Finally, the equilibration protocol provides a system ready for production MD simulations within a few days on broadly accessible GPU resources. The pipeline is applied to study the Zika virus as a test case for large biomolecular systems. The VLP stabilizes at approximately 0.5 µs of MD simulation, reproducing correlations greater than 0.90 against experimental density maps from cryo-electron microscopy. Detailed structural analysis of the protein envelope also shows very good agreement in root-mean-square deviations and B-factors with the experimental data. The level of details attained shows for the first time a possible role for anionic phospholipids in stabilizing the envelope. Combining an efficient and reliable setup procedure with an accurate coarse-grained force field provides a valuable pipeline for simulating arbitrary viral systems or subcellular compartments, paving the way toward whole-cell simulations.

Nitro-fatty acids as activators of hSIRT6 deacetylase activity.

Sirtuin 6, SIRT6, is critical for both glucose and lipid homeostasis and is involved in maintaining genomic stability under conditions of oxidative DNA damage such as those observed in age-related diseases. There is an intense search for modulators of SIRT6 activity, however, not many specific activators have been reported. Long acyl-chain fatty acids have been shown to increase the weak in vitro deacetylase activity of SIRT6 but this effect is modest at best. Herein we report that electrophilic nitro-fatty acids (nitro-oleic acid and nitro-conjugated linoleic acid) potently activate SIRT6. Binding of the nitro-fatty acid to the hydrophobic crevice of the SIRT6 active site exerted a moderate activation (2-fold at 20 µm), similar to that previously reported for non-nitrated fatty acids. However, covalent Michael adduct formation with Cys-18, a residue present at the N terminus of SIRT6 but absent from other isoforms, induced a conformational change that resulted in a much stronger activation (40-fold at 20 µm). Molecular modeling of the resulting Michael adduct suggested stabilization of the co-substrate and acyl-binding loops as a possible additional mechanism of SIRT6 activation by the nitro-fatty acid. Importantly, treatment of cells with nitro-oleic acid promoted H3K9 deacetylation, whereas oleic acid had no effect. Altogether, our results show that nitrated fatty acids can be considered a valuable tool for specific SIRT6 activation, and that SIRT6 should be considered as a molecular target for in vivo actions of these anti-inflammatory nitro-lipids.

Split the Charge Difference in Two! A Rule of Thumb for Adding Proper Amounts of Ions in MD Simulations.

Despite the relevance of properly setting ionic concentrations in Molecular Dynamics (MD) simulations, methods or practical rules to set ionic strength are scarce and rarely documented. Based on a recently proposed thermodynamics method we provide an accurate rule of thumb to define the electrolytic content in simulation boxes. Extending the use of good practices in setting up MD systems is promptly needed to ensure reproducibility and consistency in molecular simulations.

Fat SIRAH: Coarse-Grained Phospholipids To Explore Membrane-Protein Dynamics.

The capability to handle highly heterogeneous molecular assemblies in a consistent manner is among the greatest challenges faced when deriving simulation parameters. This is particularly the case for coarse-grained (CG) simulations in which chemical functional groups are lumped into effective interaction centers for which transferability between different chemical environments is not guaranteed. Here, we introduce the parametrization of a set of CG phospholipids compatible with the latest version of the SIRAH force field for proteins. The newly introduced lipid species include different acylic chain lengths and partial unsaturation, as well as polar and acidic head groups that show a very good reproduction of structural membrane determinants, such as areas per lipid, thickness, order parameter, etc., and their dependence with temperature. Simulation of membrane proteins showed unprecedented accuracy in the unbiased description of the thickness-dependent membrane-protein orientation in systems where this information is experimentally available (namely, the SarcoEndoplasmic Reticulum Calcium-SERCA-pump and its regulator Phospholamban). The interactions that lead to this faithful reproduction can be traced down to the single amino acid-lipid interaction level and show full agreement with biochemical data present in the literature. Finally, the present parametrization is implemented in the GROMACS and AMBER simulation packages facilitating its use by a wide portion of the biocomputing community.

From quantum to subcellular scales: multi-scale simulation approaches and the SIRAH force field.

Modern molecular and cellular biology profits from astonishing resolution structural methods, currently even reaching the whole cell level. This is encompassed by the development of computational methods providing a deep view into the structure and dynamics of molecular processes happening at very different scales in time and space. Linking such scales is of paramount importance when aiming at far-reaching biological questions. Computational methods at the interface between classical and coarse-grained resolutions are gaining momentum with several research groups dedicating important efforts to their development and tuning. An overview of such methods is addressed herein, with special emphasis on the SIRAH force field for coarse-grained and multi-scale simulations. Moreover, we provide proof of concept calculations on the implementation of a multi-scale simulation scheme including quantum calculations on a classical fine-grained/coarse-grained representation of double-stranded DNA. This opens the possibility to include the effect of large conformational fluctuations in chromatin segments on, for instance, the reactivity of particular base pairs within the same simulation framework.

The SIRAH 2.0 Force Field: Altius, Fortius, Citius.

A new version of the coarse-grained (CG) SIRAH force field for proteins has been developed. Modifications to bonded and non-bonded interactions on the existing molecular topologies significantly ameliorate the structural description and flexibility of a non-redundant set of proteins. The SIRAH 2.0 force field has also been ported to the popular simulation package AMBER, which along with the former implementation in GROMACS expands significantly the potential range of users and performance of this CG force field on CPU/GPU codes. As a non-trivial example of its application, we undertook the structural and dynamical analysis of the most abundant and conserved calcium-binding protein, calmodulin (CaM). CaM is composed of two calcium-binding motifs called EF-hands, which in the presence of calcium specifically recognize a cognate peptide by embracing it. CG simulations of CaM bound to four calcium ions in the presence or absence of a binding peptide (holo and apo forms, respectively) resulted in good and stable ion coordination. The simulation of the holo form starting from an experimental structure sampled near-native conformations, retrieving quasi-atomistic precision. Removing the binding peptide enabled the EF-hands to perform large reciprocal movements, comparable to those observed in NMR structures. On the other hand, the isolated peptide starting from the helical conformation experienced spontaneous unfolding, in agreement with previous experimental data. However, repositioning the peptide in the neighborhood of one EF-hand not only prevented the peptide from unfolding but also drove CaM to a fully bound conformation, with both EF-hands embracing the cognate peptide, resembling the experimental holo structure. Therefore, SIRAH 2.0 shows the capacity to handle a number of structurally and dynamically challenging situations, including metal ion coordination, unbiased conformational sampling, and specific protein-peptide recognition.

FRET biosensor uncovers cAMP nano-domains at ß-adrenergic targets that dictate precise tuning of cardiac contractility.

Compartmentalized cAMP/PKA signalling is now recognized as important for physiology and pathophysiology, yet a detailed understanding of the properties, regulation and function of local cAMP/PKA signals is lacking. Here we present a fluorescence resonance energy transfer (FRET)-based sensor, CUTie, which detects compartmentalized cAMP with unprecedented accuracy. CUTie, targeted to specific multiprotein complexes at discrete plasmalemmal, sarcoplasmic reticular and myofilament sites, reveals differential kinetics and amplitudes of localized cAMP signals. This nanoscopic heterogeneity of cAMP signals is necessary to optimize cardiac contractility upon adrenergic activation. At low adrenergic levels, and those mimicking heart failure, differential local cAMP responses are exacerbated, with near abolition of cAMP signalling at certain locations. This work provides tools and fundamental mechanistic insights into subcellular adrenergic signalling in normal and pathological cardiac function.

Regulation of signaling directionality revealed by 3D snapshots of a kinase:regulator complex in action.

Two-component systems (TCS) are protein machineries that enable cells to respond to input signals. Histidine kinases (HK) are the sensory component, transferring information toward downstream response regulators (RR). HKs transfer phosphoryl groups to their specific RRs, but also dephosphorylate them, overall ensuring proper signaling. The mechanisms by which HKs discriminate between such disparate directions, are yet unknown. We now disclose crystal structures of the HK:RR complex DesK:DesR from Bacillus subtilis, comprising snapshots of the phosphotransfer and the dephosphorylation reactions. The HK dictates the reactional outcome through conformational rearrangements that include the reactive histidine. The phosphotransfer center is asymmetric, poised for dissociative nucleophilic substitution. The structural bases of HK phosphatase/phosphotransferase control are uncovered, and the unexpected discovery of a dissociative reactional center, sheds light on the evolution of TCS phosphotransfer reversibility. Our findings should be applicable to a broad range of signaling systems and instrumental in synthetic TCS rewiring.

SIRAH tools: mapping, backmapping and visualization of coarse-grained models.

UNLABELLED: Coarse-grained (CG) models reduce the cost of molecular dynamics simulations keeping the essence of molecular interactions. Still, the diversity of CG representations (sizes, connectivity, naming, etc.) hampers the handling and visualization of such models. SIRAH Tools comprises a set of utilities to convert all-atoms coordinates to arbitrary residue-based CG schemes, write GROMACS' topological information at any resolution into PSF format and a VMD plugin to visualize, analyze and retrieve pseudo-atomistic information from CG trajectories performed with the SIRAH force field. These tools facilitate the use of intricate CG force fields outside the small developer's community. AVAILABILITY AND IMPLEMENTATION: Different utilities of SIRAH Tools are written in Perl, Tcl, or R. Documentation and source codes are freely distributed at http://www.sirahff.com CONTACT: : mmachado@pasteur.edu.uy or spantano@pasteur.edu.uy.

Exploring LacI-DNA dynamics by multiscale simulations using the SIRAH force field.

The lac repressor protein (LacI) together with its target regulatory sequence are a common model for studying DNA looping and its implications on transcriptional control in bacteria. Owing to the molecular size of this system, standard all-atom (AA) simulations are prohibitive for achieving relevant biological time scales. As an alternative, multiscale models, which combine AA descriptions at particular regions with coarse-grained (CG) representations of the remaining components, were used to address this computational challenge while preserving the relevant details of the system. In this work, we implement a new multiscale approach based on the SIRAH force field to gain deeper insights into the dynamics of the LacI-DNA system. Our methodology allows for a dual resolution treatment of the solute and solvent, explicitly representing the protein, DNA, and solvent environment without compromising the AA region. Starting from the P1 loop configuration in an undertwisted conformation, we were able to observe the transition to the more stable overtwisted state. Additionally, a detailed characterization of the conformational space sampled by the DNA loop was done. In agreement with experimental and theoretical evidence, we observed the transient formation of kinks at the loop, which were stabilized by the presence of counterions at the minor groove. We also show that the loop's intrinsic flexibility can account for reported FRET measurements and bent conformations required to bind the CAP transcription factor.

Isoform-specific determinants in the HP1 binding to histone 3: insights from molecular simulations.

Despite the significant improvements in anti HIV-1 treatment, AIDS remains a lifelong disease due to the impossibility to eradicate the viral reservoir established upon integration of the viral genome. Controlling the epigenetic block imposed by the host cell machinery to the viral transcription may represent a therapeutic alternative to purge the viral reservoir, offering a way to eradicate the infection. Heterochromatin protein 1 (HP1) has been reported to actively participate in the silencing of HIV-1 integrated genome by binding to histone 3 (H3) tail. This interaction is mediated by the Chromodomain of HP1. Nevertheless, the structural features that determine its binding to H3 tail upon post-transductional modifications, such as methylation and phosphorylation as well as isoform-specific effects have not yet been described. We have undertaken the systematic simulation of the Chromodomains of the isoforms beta and gamma of HP1 in complex with the H3 tail methylated at Lys9 in presence/absence of phosphorylation at Ser10. Our results pinpoint isoform-specific electrostatic interactions as important determinants for the stability of the complexes. Characterization of intermolecular contacts between HP1 variants and H3 furnishes new insights on isoform-specific recognition and the effect of phosphorylation.

A Coarse Grained Model for Atomic-Detailed DNA Simulations with Explicit Electrostatics.

Coarse-grain (CG) techniques allow considerable extension of the accessible size and time scales in simulations of biological systems. Although many CG representations are available for the most common biomacromolecules, very few have been reported for nucleic acids. Here, we present a CG model for molecular dynamics simulations of DNA on the multi-microsecond time scale. Our model maps the complexity of each nucleotide onto six effective superatoms keeping the "chemical sense" of specific Watson-Crick recognition. Molecular interactions are evaluated using a classical Hamiltonian with explicit electrostatics calculated under the framework of the generalized Born approach. This CG representation is able to accurately reproduce experimental structures, breathing dynamics, and conformational transitions from the A to the B form in double helical fragments. The model achieves a good qualitative reproduction of temperature-driven melting and its dependence on size, ionic strength, and sequence specificity. Reconstruction of atomistic models from CG trajectories give remarkable agreement with structural, dynamic, and energetic features obtained from fully atomistic simulation, opening the possibility to acquire nearly atomic detail data from CG trajectories.